Dispersion properties of transverse waves propagating in the electrically polarized BEC 
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Further development of the method of the quantum hydrodynamics in application to a Bose-Einstein 
condensate is presented. To consider evolution of polarization direction along with particles movement we 
have developed the set of quantum hydrodynamics equations, which includes equation of the polarization 
evolution and the polarization current evolution along with the continuity equation and the Euler equation 
(momentum balance equation). Dispersion properties of transverse waves including the electromagnetic 
waves propagating through the BEC are considered. For this purpose we consider whole system of the 
Maxwell's equations for field dynamic description. This approximation gives us possibility to consider 
electromagnetic waves along with matter waves. We show that in the electrically polarized BEC we have 
a splitting of the electromagnetic waves on two branches. We have four solution, two for the electro- 
magnetic waves and two for the matter waves, the last two are the concentration-polarization waves. In 
this approximation we find anisotropy in dispersion dependencies for the all four waves. Whereas in the 
longitudinal approximation we have two matter waves only which show no anisotropy. 



I. INTRODUCTION 

At studying of an electrically polarized Bose-Einstein 
condensate (BEC) the generalization of the Gross- 
Pitaevskii (GP) equation [Qj]- ]5|] and some it's further gen- 
eralization |@], 0], see also review papers JUl- lUIll . es- 
pecially we want to admit recent review written by M. A. 
Baranov et. al. [ .12]]. This nonlinear Schrodinger equation 
defines the complex scalar wave function ip(r,t), which, 
as in the case of the non-polarized BEC, describes evo- 
lution of concentration n(r, t) = ^*(r,t)ip(r,t). When 
direction of polarization of each particle d is not changing 
during the considering process and all of them have same 
direction we can write formula for density of polarization 
P = ip*(r, t)dip(r, t). Thus, polarization density changes 
due to particles movement, and, consequently, changing of 
concentration. The waves of polarization can easily prop- 
agate even in the ferroelectrics as well as spin waves in 
ferromagnetics. All the more we can expect existence of 
them in polarized BEC. To consider evolution of polar- 
ization direction along with particles movement we have 
developed the set of quantum hydrodynamics equations, 
which includes equation of the polarization evolution and 
the polarization current evolution along with the continuity 
equation and the Euler equation (the momentum balance 
equation) Q- 0]. 



It was shown in Ref. 111311 that Hamiltonian of the elec- 
tric dipole interaction correct form which corresponds to 
the Maxwell's equation is 
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Using well-known identity 
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we can see that the Hamiltonian {T} differs from usually 
used one 
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Necessity to consider the Hamiltonian of electric dipole in- 
teraction in the form (fl} caused by the fact that it must 
accord to the Maxwell's equation. This connection exists 
since Maxwell's equations describe electric field and it's 
connection with the sources. In our case source is the den- 
sity of electric polarization. Action of the electric field on 
the density of polarization come in to equations of motion 
via the force field. As the result it describes interaction of 
polarization (electric dipole moments). 

It has been shown in Ref.s lfl7ll - Ill9ll that using of the 
non-linear Schrodinger equation which generalized the GP 
equation for the system of the electrically polarized BEC 
and using of the formula OJ for description of the dipole - 
dipole interaction lead to the appearance of the contribu- 
tion of the equilibrium polarization in the dispersion de- 
pendence of the Bogoliubov's mode, this contribution leads 
to the anisotropy of the dispersion dependence. In contrast 
with it, in our previous papers lfl3ll - iflfjh . we have shown 
that at account of the Hamiltonian £Q and consideration of 
evolution of electric dipole directions we have got disper- 
sion dependencies for two waves, that means what account 
of polarization evolution leads to appearing of new wave 
in the BEC. Moreover, we have found that these dispersion 
dependencies are isotropic, and contain no dependence on 
direction between direction of the wave propagation and 
direction of the equilibrium polarization. 

In this paper we suggest more general model in com- 
parison with our recent papers. We consider whole system 
of the Maxwell's equations instead of the pair quasi-static 
equations describing the electric field considered earlier. 
This approximation gives us possibility to consider elec- 
tromagnetic waves along with matter waves. Below we 
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will show that in the electrically polarized BEC we have 
a splitting of the electromagnetic waves on two branches. 
So, we have four solution, two for the electromagnetic 
waves and two for the matter waves, the last two are waves 
of concentration-polarization. It is important and interest- 
ingly to admit that in this approximation we find anisotropy 
in the dispersion dependencies for all four waves. 

This paper is organized as follows. In Sec. II we present 
the quantum hydrodynamics model for description of the 
electrically polarized BEC. In Sec. Ill we present disper- 
sion equation and describe the method to get it. In. Sec. 
IV we consider important consequence of the obtained dis- 
persion equation. In Sec. V brief discussion of obtained 
results is presented. 



H. BASIC EQUATIONS 

In our previous paper lfl3ll we have presented the de- 
tailed derivation of the QHD equation for the electrically 
polarized BEC. These equations were directly derived from 
the microscopic many-particle Schrodinger equation. Let 
us briefly present the set of the quantum hydrodynamics 
equations. 

The first equation of the QHD equations system is the 
continuity equation 

d t n + d a {nv a ) = 0. (4) 

The momentum balance equation for the polarized BEC 
has the form 
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are the numerical coefficients. In equation (O we define a 
parameter T 2 as (0. This definition differs from the one in 
the work 12011 . Terms proportional to h 2 appear as a result 
of usage of the quantum kinematics. The first two members 
at the right side of the equation (Q are first terms of expan- 
sion of the quantum stress tensor. They occur because of 
taking into account of the SRI potential The interac- 
tion potential determinates the macroscopic interaction 
constants T and T 2 . The last two members of the equa- 
tion (0 describe force fields that affect the dipole moment 



in a unit of volume as the effect of the external electrical 
field and the field produced by other dipoles, respectively. 
The last member is written using the self-consistent field 
approximation iflol . 12111 . p a/3 (r,t) is the tensor of the 
kinetic pressure, which depends on particle thermal veloc- 
ities and makes no contribution into the BEC dynamics at 
temperatures near zero. 

The first order by the interaction radius constant for di- 
lute gases has the form 
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where a is the scattering length SL 
We have also the field equations 

VE(r,t) = -4vrVP(r,t), 

and 

V x E(r,t) = 0. 

In the case particles does not contain the dipole moment, 
the continuity equation and the momentum balance equa- 
tion form a closed system of equations. The new physical 
value emerges, when the dipole moment is taken into ac- 
count in the momentum balance equation, it is the polariza- 
tion vector field P a (r, t). This causes system of equations 
to become incomplete. 

Next equation we need for investigation of the collective 
excitation dispersion is the equation of polarization evolu- 
tion 



0. 
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i? a/3 (r, t) is the current of polarization. 

The equation (fTTT l does not contain information about 
the effect of the interaction on the polarization evolution. 
The evolution equation of R a ^(r, t) can be constructed by 
analogy with the above derived equations. Using the self- 
consistent field approximation of the dipole-dipole inter- 
action we obtain the equation for the polarization current 
R a0 (r,t) evolution 
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Here r a/37 (r, t) presents the contribution of the thermal 
movement of the polarized particles into the dynamics of 
i? a/3 (r, t). As we deal with the BEC below, the contribu- 
tion of r a ^ 7 (r, t) may be neglected. The last term of the 
formula (fT2l) includes both external electrical field and the 
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self-consistent field that particle dipoles create. This term 
contain numerical constant a. 

Equations Q and ( fTOt emerge in the considered non- 
relativistic limit and allow to study longitudinal waves 
only. When we use term longitudinal we suggest that di- 
rection of the electric field perturbation is parallel to the 
direction of the wave propagation. This case has been con- 
sidered in Ref.s lfl3ll - lll5ll It is well-known that equa- 
tions © and ([Tol l are the part of the set of the Maxwell's 
equations. Therefore, if we want to consider the transverse 
waves we should use the whole set of the Maxwell's equa- 
tions, which has well-known form 



E°> = E% + 5E a , B a = + 8B C 



pa = pa + gp* ficfi = + gfiat 



(17) 



VB(r,t) = 0, 



(13) 



Substituting these relations into system of equations ©, 
©, (fTTT l. (fT2l) and (fT3l- (fT6l) and neglecting nonlinear 
terms, we obtain a system of the linear homogeneous equa- 
tions in partial derivatives with constant coefficients. Pass- 
ing to the following representation for small perturbations 
Sf 

Sf = f(u,k)exp(—iujt + ikr) 



and 



VE(r,t) = -4vrVP(r,t), 
VxE(r,i) = --d t B(r,t), 



(14) 
(15) 
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yields the homogeneous system of algebraic equations. 
The electric field strength is assumed to have a nonzero 
value. Expressing all the quantities entering the system 
of equations in terms of the electric field, we come to the 
equation 



where 
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III. ELEMENTARY EXCITATIONS IN THE POLARIZED 
BEC: DISPERSION EQUATION 



k 2 ]S ap + k a k p 



We can analyze the linear dynamics of collective excita- 
tions in the polarized BEC using the QHD equations (O, 
©, £□}, £E2]l and £[3])-£[6]). Let's assume the system is 
placed in an external electrical field E = E e z . The val- 
ues of concentration n and polarization P = kEq for 
the system in an equilibrium state are constant and uniform 
and its velocity field v a (r, t), magnetic field B a (r, t) and 
tensor R a/3 (r, t) values are zero. 

We consider the small perturbation of equilibrium state 

n = n + 6n, v a = 4- 8v a , 
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where k 2 = kl + k 2 + k 2 . 

- L y z 

We suppose that the external electric field parallel to the 
z-axis. Thereby, the equilibrium polarization P also par- 
allel to the z-axis P || e z . In this case tensor ( fT9l l become 
more simple 
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Dispersion equation is detA = 0. Using formula (j20|) 
we have evident form of dispersion equation 

(u 2 - k 2 c 2 ) 2 

+ (3(oj)(lu 4 - (1 + cos 2 9)oj 2 k 2 c 2 + cos 2 9k 4 c 4 ) = 0. 

(21) 



IV. DISPERSION DEPENDENCIES 

In the low frequency limit (oj <C kc) from equation (l2ll 
we find 

1 + cos 2 6P(u) = 0, (22) 
at cos 2 9 = 1 we have 

1 + = 0. (23) 



Equation (l23l was obtained in Ref.s II 1 311 . II 1411 . II 1 511 and 
was in detail studied in Ref. Ill 311 . 
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If P Q = we find ui = kc. 
Equation d22l has following solution 
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In the case of wave propagation parallel to the equilibrium 
polarization cos 9 = 1 solution we get solution obtained 
in Ref. s lfl3ll . [14], Il5ll . and at great length considered in 
Ref. [13]. At increasing of angle 9 from to tt/2 con- 
tribution of equilibrium polarization become smaller and 
vanish at 9 = tt/2. Obtaining of an anizotropic dispersion 
dependence emerges because of consideration of the whole 
set of the Maxwell's equations (fl"3l-(fT6l). It gives us possi- 
bility to account transverse electric field in the wave, along 
with the longitudinal field, the last one was considered in 
our previous papers lfl3ll . ifbUl . lfl5ll . 

Under the condition 9 = or 9 = ir equation (1211 as- 
sumes the form 



(w 2 - A: 2 c 2 ) 2 (l + P(u)) = 0. 



(25) 



From this equation we find usual the dispersion depen- 
dence of the light to = kc, and equation 1 + /3(cj) = 
describes the two matter waves, which appear instead of 
the Bogoliubov's mode in the electrically polarized BEC, 
see Ref.s HI, Q, JH. 

In the opposite case, when 9 = tt/2 equation (I2H sim- 
plifies to 



(to 2 



k 2 c 2 ){co 2 - k 2 c 2 + P(to)to 2 ) = 0. 



(26) 



In this case we have two independent equation. One of 
them to 2 — k 2 c 2 = describes changeless dispersion of 
the light. The second equation 



w 2 (l + P(u)) - k 2 c 2 







(27) 



is the equation of the third degree of to 2 . So, it contains 
three wave branches. We can admit that it contains both 
modified by the BEC dispersion of the light and the two 
matter wave. This equation has no low frequencies limita- 
tion as d22l . thus we have here high frequencies "tail" in 
the dispersion of the matter waves. 

Let's consider influence of the BEC on the light as a 
small correction. Therefore we put to = kc in f3(co) 
of equation (l27l l. and accounting that kc >> hk 2 /m, 



^Tn /mk, y/T 2 n /mk 2 , in the result we find 



to ~ kc/Jl + /3(kc), 



(28) 



where /3(kc) = — 2ttTPq /(m 2 n Q kc 3 ). Expanding for- 
mula (l28l) in a series by (3(kc) we get 



to 
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We can see that the polarized BEC caused splitting of the 
light on two waves at light propagation perpendicular to 
the direction of the equilibrium polarization. Magnitude of 
splitting M s is equal to 



ttTP 2 

m 2 n c 2 
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V. CONCLUSION 

We have considered transverse waves in the electrically 
polarized BEC. To do it we have used whole set of the 
Maxwell's equation for description of electromagnetic field 
caused by the dynamic of electric dipoles of the medium. 
In comparison to the case of longitudinal waves considered 
in our previous papers we have found two advantages. 

The first of them is the possibility of considering of light 
propagating through the electrically polarized BEC and 
studying of it's properties. Thus, we have considered dis- 
persion of the light propagating through the medium and 
we have found the splitting of the light on two waves in 
the case when the light propagate at the angle to the direc- 
tion of the equilibrium polarization. We have calculated 
the magnitude of frequencies splitting for the case of the 
light propagation perpendicular to the direction of equilib- 
rium polarization. There is no splitting in the case of the 
light propagation along the direction of the equilibrium po- 
larization, there is also no contribution of the medium in 
the light dispersion dependence in this case. 

The second advantage of the considered approximation 
is the appearing of the anisotropy in the dispersion depen- 
dencies of the matter waves. As in the case of longitudinal 
waves we have got two waves which exist instead of the 
Bogoliubov's mode existing in the unpolarized BEC. How- 
ever, we did not find, in our previous papers, dependence 
of the frequencies of the collective excitation on the direc- 
tion of wave propagation for the longitudinal waves, which 
might be expected due to the anisotropy of the Hamil- 
tonian of dipole-dipole interaction. We have found ana- 
lytical solution for the plane matter waves propagating in 
the three dimensional BEC and shown the influence of the 
anisotropy on the value of the frequencies. 

In the result, we have developed generalization of the 
quantum hydrodynamic equations for the electrically po- 
larized BEC including transverse wave propagation. Using 
it we have shown that in the electrically polarized BEC sit- 
uated in the external electric field there are four wave, two 
of them are high frequencies waves, and them associated 
with the light. Two more wave are the anisotropic matter 
waves. 
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